Running title: Compact sources in M 82 



Could the Compact Radio Sources in M 82 be Cluster Wind 

Driven Bubbles? 

y^Q ] E. R. Seaquist and M. Stankovic 

Department of Astronomy & Astrophysics, University of Toronto, 50 St. George Street, 

Toronto, Ontario, M5S 3H4, Canada 



O 
<D 

Q- seaquist@astro.utoronto.ca 

ABSTRACT 

> \ 

The compact non-thermal sources in M82 and other starburst galaxies are 
generally thought to be supernova remnants (SNRs). We consider an alterna- 
tive hypothesis that most are wind driven bubbles (WDBs) associated with very 
young super star clusters (SSCs). In this scenario, the synchrotron emitting par- 
ticles are produced at the site of the shock transition between the cluster wind 
and the hot bubble gas. The particles radiate in the strong magnetic field pro- 
duced in the expanding shell of shocked ambient interstellar gas. One of the 
+3 . motivations for this hypothesis is the lack of observed time variability in most 

of the sources, implying ages greater than expected for SNRs, but comfortably 
within the range for WDBs. In addition, as SNRs, these sources are not effective 
in driving the starburst mass outflow associated with the nuclear region of M82, 
thus requiring a separate mechanism for coupling SN energy to this outflow. The 
WDB hypothesis is found to be feasible for underlying clusters in the mass range 
~ 2 x lO 4 ^ 1 ) M and ambient gas densities in the range ~ 3 x 1 3(±D cm . 
The ages of the bubbles are between several xlO 3 and several xlO 4 years. Since 
the SNR picture cannont be ruled out, we provide suggestions for specific ob- 
servational tests which could confirm or rule out the WDB hypothesis. Finally, 
we discuss the WDB hypothesis in the context of broader phenomena in M82, 
such as the rate of star formation and starburst outflows, and the possible inter- 
pretation of supershells in M82 as the products of multiple supernovae in young 
SSCs. 



Subject headings: galaxies: individual (M82) — stars: winds, outflows — super- 
nova remnants 
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INTRODUCTION 



The nearby starburst galaxy M82 is host to a population of at least 60 compact (di- 
ameter a few pc) radio sources observable at centimeter wavelengths. Similar populations 
may be found in a number of other starburst systems including the ultra-luminous IR galaxy 
Arp220. For the approximately 60 sources in M82, about half are identified (by their spec- 
tral indices) as non-thermal or synchrotron emitting sources, and are generally accepted to 
be the remnants of supernovae associated with individual stars. About a quarter are identi 



fiedas thermal sources or compact HII regions ([McDonald et alJl2002j ; iRodriguez-Rico et al. 



20041 ) and the rest are not yet identified. The association of the compact non-thermal ra- 
dio sources (CNRs) with supernova remnants (SNRs) is based on the similarity of their 
luminosities, morphology and nonthermal spectra to those of SNRs in our Galaxy. 

The purpose of this paper is to explore an alternative interpretation that many, and 
possibly most, of the CNRs are generated by wind driven bubbles (WDBs) associated with 
very young super star clusters (SSCs). In this picture, the strong stellar wind from a young 
cluster energizes a hot expandi ng bubble which swe e ps up interstellar ma terial to form a thin 



and dense radiative shell (e.g. IWeaver et al.lll977l ; iKoo fc McKed 119921 ) . In the suggestion 



presented in this paper, the non-thermal emission is associated with relativistic electrons 
accelerated at the adiabatic shock responsible for heating the bubble. The non-thermal 
emission occurs as these electrons interact with the interstellar magnetic field compressed in 
the radiative shell comprising the shocked interstellar medium (ISM). Thus the morphology 
of this WDB would indeed resemble that of an SNR. 

But why consider this alternative if the SNR hypothesis works? The motivation is 
driven by some troubling, if non-fatal, issues with the SNR hypothes is. First, a majority of 
the sources show no variations in flux density. iKronberg et al.l (120001 ) have studied the time 
variability of a sample of 24 sources, and have provided a stringent upper limit of 0.1% per 
year for the degree of variability in a sub-sample of 18 over 12 years. This suggests a lower 
limit on the age of at least 1000 years. If, on the other hand, SNRs with diameters up to 4pc 
(comparable to those in M 82) are i n the free expansion stage with velocities ^ 5000 km s -1 , as 



suggested by iMuxlow et al.l (119941 ) . the corresponding age would be at most a few hundred 
years. However the SNRs in M82 could be decelerated by dense gas, which would make 
them significantly older. In this case the inferred ages greater than 1000 years correspond 
to expansion velocities <1000kms _1 , and lower than this at later stages. It is difficult to 
see how such supernovae (SNe) could drive the starburst nuclear outflow at several thousand 
kms" 1 . Secondly, a significant fraction (perhaps most?) of the massive stellar p opulation 
is formed in SSCs, as suggested by some authors (e.g. iLipscy fc Plavchanl 120041 ) . Within 
such clusters, SNe would be frequent and collectively produce a nearly steady wind rather 
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than individually identifiable SNRs (IMcCrav fc Kafato Jl987l : Mac Low & McCravill988[ ) . The 



result would be shells that look very large and diffuse compared with those for individual 
SNRs. We return to this point later to discuss evidence for such structures in M82. 



Some of these difficulties were recently addressed by I Chevalier fc Franssonl (120011 ). To 
account for the lack of variability, they suggested that most of the SNRs have recently en- 
tered the radiative phase as they expand into an ISM of density 10 3 H atoms cm -3 . The 
corresponding expansion velocities are ~500kms _1 . To drive the starburst wind, they pos- 
tulate the existence of another population of SNRs expanding into a lower density medium, 
permitting mechanical energy to be transferred to the starburst wind. Thus the SNR model 
is certainly viable, but in light of the concerns which have spawned such discussion, the 
WDB hypothesis should be considered. 

The WDB hypothesis predicts ages of at least several thousand years for the compact 
sources, much greater than for SNRs. Figure 1 shows a comparison between the evolution 
of an SNR (both adiabatic and radiative), and a WDB, both expanding into a medium 
of density 10 3 cm~ 3 . The WDB curves correspond to SSC masses of 10 4 and 1O 5 M , not 
atypical for M82. The WDB evolution represented are described by equations (1) and (2) 
in § 2, where the evolution is discussed in more detail. The point to be made in Figure 1 
is simply that for the typical radii of the non-thermal shells (0.5-2pc), WDBs are much 
older than SNRs. In the WDB picture, the modest shells produced by the stellar winds are 
replaced after several million years by supershells driven by the collective winds of many SNe 
within the cluster, as noted earlier. It is the latter winds that are collectively responsible for 
energizing the starburst outflow. 

We concern ourselves in this paper with an evaluation of the WDB picture as the 
explanation for the CNRs, using existing data in the literature. It is not the intention to 
propose that the WDB hypothesis should now replace the conventional SNR picture, but 
rather to determine its feasibility by subjecting it to a number of tests. The general approach 
taken in the paper is to assume that the CNRs have a WDB origin, and then outline and 
discuss the consequences. 

In § 2, we outline the model. Section 3 is devoted to describing tests to determine the 
feasibility of the model and ranges in parameter space for the viability of the model. The 
results are presented in § 4. In § 5, we present inferences from these results concerning the 
properties of the shells, and in § 6 we discuss the relativistic particle acceleration efficiencies 
with a comparison with Galactic SNRs. Section 7 addresses the issue of the source sizes 
compared to the radii of the hypothesized underlying clusters, and § 8 assesses the likelihood 
of detecting these clusters and their effects in various wavebands. Finally we present in § 9 the 
conclusions, placing the WDB hypothesis in a larger context of other starburst phenomena 
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Fig. 1. — Plots of radius vs. time for SNR and WDB shells. The two curves representing 
SNRs show adiabatic and radiative behavior, and the two curves representing WDBs are 
based on winds for cluster masses of 10 4 and 10 5 M Q . The ISM density is 10 3 cm -3 in all 
cases. 

in M82, and outlining stringent tests for the WDB picture. 



2. OUTLINE OF THE MODEL 



The picture underlying our hypothesis is based on the dynamical evolution of a hot 
WDB expand i ng in to a uniform surrou n ding ISM, as developed i n a s eries of papers by 
Weaver et all Jl977h . iMcCrav fc Kafatosl Jl987h and iKoofcMcKeel Jl992h . In this case the 
wind is the combined winds of all massive stars in a very young SSC, and the ISM is either 
the ambient density of the molecular cloud containing the SSC, the ambient ISM, or possibly 
the remnant medium produced by bi-polar flows in the proto-stellar stage of the cluster stars. 
For simplicity we overlook possible departures in the bubble geometry from that described 
in the above papers attributable to the finite size of the cluster. The expanding bubble 
drives a shock into the surrounding ISM forming a thin radiative shell containing swept up 
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interstellar gas and magnetic field. The radius r varies according to the similarity solution, 



r = 0.0170 (L 38 /n ) 1/5 t 3/5 pc 



where t is the time (yr), L38 is the wind power in units of 10 38 ergs _1 , and no is the 
uniform ISM density in H atoms cm -3 . This equation is analogous to the Sedov adiabatic 
blast wave similarity solution relevant to an adiabatically expanding SNR, given by 



r = 0.310 (E 51 /n ) 1/5 t 2/5 pc 



(2) 



where E 51 is the mechanical energy associated with the blast ejecta in units of 10 51 erg. 

Equations (1) and (2) are represented in Figure 1, except for the addition of a curve 
for the radiative phase of SNR expansion which simply conserves momentum. For the SNR 
behavior we have used E51 = 1 and n = 10 3 cm~ 3 . The radiative phas e is assumed to begin 
at r = 1 pc, appropriate for an ISM density of 10 3 cm -3 (see for e.g. I Chevalier fc Fransson 



200ll ) . The wind power level s in the figure are based on models from Starburst99 (hereafter 
SB99) (ILeitherer et al.lll999l ). These levels are presented in §3. 



The interior of the wind bubble comprises a rarified hot (T ~ 10 6 K) plasma, the conse- 
quence of a shock transition between the hypersonic stellar wind and the gas in the bubble 
cavity (see also section 6.3). The hot bubble is surrounded by a contact discontinuity sepa- 
rating the bubble from the shocked ISM. In this model the bubble shock is the source of the 
relativistic electrons causing the observed synchrotron radiation. The relativistic particles 
are assumed to propagate toward the contact discontinuity given by equation (1), where 
they diffuse into the region of the radiative dense shell and emit synchrotron radiation in the 
compressed magnetic field of the ISM. In addition, these electrons may be further accelerated 
at the interstellar shock. 

The mechanism for accelerating the relativistic particles is diffusive shock accelera- 
tion, or first order Fermi mechanism, involving particles trapped between collapsing mirrors 
formed by the bubble shock and upst ream magnetic irregularities (see for e.g. reviews by 
Drurylll983l ; iBlandford fc Eichlerl 119871 ) . The mechanism is identical to that thought to pro- 
duce the synchrotron emitting electrons in SNRs and many other relevant astrophysical sce- 
nari os. We note in part icular the shock induced synchrotron e mission in Wolf-Rayet binaries 



Pittard et al. hood ) and shocks embedded in O star winds ( iVan Loo. Runacres fc Blommd 



20061 ). One of the scenarios described in the aforementioned reviews is particle acceleration 
at the termination shock of a stellar wind, which is essentially the nature of the bubble 
shock. The treatment of the acceleration process is beyond the scope of this paper. Instead, 
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the production of non-thermal electrons involving shock speeds of several thousand kms -1 
is taken as assured, since it occurs for example in SNRs. The issue then focuses on the 
acceleration efficiency, which is discussed in §6. 

The model implicitly assumes that the accelerated particles reach the outer emitting 
region in a time short compared with the age of the shell. Though the particles will tend to 
stream outward from the acceleration zone at relativistic speeds, they would undoubtedly 
be scattered by magnetic irregularities in the hot bubble. The material in the bubble is 
composed primarily of gas (with embedded fields) evaporated from the cool outer shell 



([Weaver et al.lll977l ). The treatment of this problem is outside the scope of this preliminary 
paper, but we note that the magnetic irregularities may take the form of Alfven waves 
produced by gas turbulenc e or generated b y the collective behavior of the streaming energetic 



particles themselves (e.g. IWentzel I Il974l ). Under these conditions, the limiting speed of 
propagation of the relativistic particles is the Alfven speed of the waves. If the energy 
density associated with these waves is comparable to the thermal energy density in the hot 
bubble, then the Alfven speed is comparable to the sound speed (and thus the speed of the 
stellar wind). Such conditions would lead to a propagation time which is short compared to 
the age of the WDB. 



3. FEASIBILITY TESTS 



The model outlined in § 2 was applied to published radio data for a sample of well 
observed CNRs in M 82 to determine for each source its feasibility and resulting constraints 
on SSC mass and ISM density (see §4). This sample is hereafter referred to as the input 
sample for the array of tests. The output sample comprises all sources from the input sample 
which pass all of the tests, yielding for each member a range of SSC mass and ISM density. 
The members of the output sample are shown in Table 1. We outline below the tests which, 
by their stringency, act to most severely test and constrain the model. Some of the tests 
invoke physical plausibility (e.g. § 3.1), whereas others refer to limits set by observations (e.g. 
§3.5). Both forms of test combine to impose limits on the cluster mass and/or the ambient 
gas density surrounding the bubble. Some further less stringent constraints, additional 
insights, and further discussion of plausibility are deferred to §5-8. The sole selection 
criterion for membership in the input sample is the requirement for observed flux densities 
at a minimum of four wavelength s, and thus w e ll defi n ed non-therma l spect ra l indic e s. Th e 
sam ple was drawn fr o m da ta in iHuang et al.l (119941 ). iMuxlow et al.l (119941 ). lAllenl (119991 ). 



and iMcDonald et al.l (120021 ). Excluding the strongly variable candidate 41.9 + 57.5, the 



input sample comprised 27 sources. Note that the input source sample contains a number 
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qf sources whose measured flux variabi l ity ( Kronberg et al.l 120001 ) and / or expansion rates 
( IMuxlow et al.ll2005l ; iFenech et al.ll2005l ; iBeswick et al.l 120061 ) are consistent with or suggest 
the SNR hypothesis. If such sources pass the WDB feasibility tests, then the tests do not 
discriminate between the models. The tests applied are outlined individually in § 3.1 through 
3.5. 



3.1. Energy conservation 

Conservation of energy requires the wind to supply sufficient mechanical energy over 
its age for the requirements of the synchrotron sources. The rel ativistic part i cle an d mag- 
netic field energies were derived from equipartition arguments in iPacholczyki (jl970l ). which 
provide minimum energy requirements. The required flux densities, spect r al indices, and 



angula r sizes f or ind ividual sources were ob t ained from iHuang et al.l (11994 ). IMuxlow et al. 



( 11994J ) , lAllenl (119991 ) , and McDonald et al.l (120021 ) . We assumed spherical geometry with 
source diameters equal to the average measured along major and minor axes. For a few 
sources with out published angular sizes we adopted an angular diameter of 0.2", following 
Alien! (Il999l ). Radio luminosities were computed over the standard frequency range 0.01 to 
100 GHz, and the distance adopted for M82 was 3.25 Mpc. The total equipartition energy 
may be expressed in the form E p = Ej (1 + k) 4//7 <£> 3//7 erg, where k and $ are respectively the 
(unknown) ratio of proton to electron energy and volume filling factor for the synchrotron 
source. Ei is then the minimum energy requirement when (1 + k) 4//7 $ 3//7 = 1. Note that the 
equipartition argument requires that the energy be divided approximately equally between 
relativistic particles and magnetic field energy. The quantity Ei was computed for each 
source, resulting in a median value for the output sample of 3.8 x 10 48 erg and a range from 
1.4 x 10 48 to 7.1 x 10 48 erg. The result for each source in the output sample is shown in Table 
1, which also shows the equipartition magnetic field Hi x [(1 + k)/0] 2//7 with the coefficient 
of Hi set to unity. 

The mass and energy outflow rates of the cluster wind are assumed to be those asso- 
ciated with the beginning of the main sequence phase. It is possible that outflows could 
be associated with massive protostars, but the outflow behavior for this phase is not well 
understood or presently computable. The wind parameters were derived using SB99 models 
assuming a cluster with an IMF with low and upper mass cutoffs of 0.1 and 100 M re- 
spectively and a break in slope at 0.5 M . Above the break, the slope is —2.30, and below 
the break, the slope is —1.30. The derived values for the wind and luminosity, with their 
dependence on SSC mass M 4 (in units of 10 4 M ) are then: 



Mass loss rate (dM/dt) = 1.86 x lO~ 5 M 4 M yr" 



(3a) 
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Wind power (L w ) = 5.61 x 10 37 M 4 ergs _1 (3b) 

Wind velocity (V w ) = 3000 km s -1 (3c) 

Bolometric luminosity (L bo i) = 3.40 x 10 40 M 4 ergs -1 (3d) 

Lyman continuum luminosity (L y ) = 4.14 x 10 50 M 4 photons s -1 (3e) 

The requirement adopted for energy conservation is that Ei/E w = f (^ 1), where E w = L w t, 



and t is the age of the wind. The age t is determined for any value of M 4 and n by 
equations (1) and (3b). The parameter f thus expresses the efficiency of conversion of wind 
energy to relativistic particle and magnetic field energy. The tests were made using f = 0.1, 
corresponding to a conversion efficiency into particle and magnetic field energy of 10% if the 
factor (l+k) 4 / 7 $ 3 / 7 = 1.0. The choice f = 0.1 is discussed and justified in §6.1. 

3.2. Radiative energy losses by relativistic electrons 

The relativistic electrons are subject to losses by synchrotron radiation and inverse 
Compton scattering. The Compton scattered photons in this case are from the cluster stars. 
The time scale or half-life of these two effects for electrons radiating at v (GHz) are given 
by 

Tsync = 0.61 H -3 / 2 i/ -1 / 2 yr (4a) 
T IC = 1.74 x 10 7 H 1/2 r 2 L 40 1 z/- 1/2 yr. (4b) 

Here H is the magnetic field strength (gauss), r is the radius of the emitting shell (pc), 
and L 40 is the bolometric luminosity of the cluster (10 40 ergs -1 ). The test requires that 
Tsync, Tic > t for v = 5 GHz to avoid significant depletion of particle energy in the observed 
region of the radio spectrum. Such depletion would produce a break in the radio spectrum 
at 5 GHz, which is not observed. The magnetic field adopted is the equipartition field that 
depends on the unknown parameters (k, $) as [(l+k)/$] 2 / 7 . This factor was set to unity in 
the computation of the quantity H used in equations (4a) and (4b). Note that the dependence 
on this unknown factor is rather weak, especially in equation (4b), where Tic ~ [(l+k) . 

We note that other types of particle losses may occur, including expansion and ionization 
losses, but these are less certain, and consideration of these is deferred to §6.3. 

3.3. Radiative losses by the hot bubble 

The initial expansion of the bubble is adiabatic, expanding according to equation (1), 
but when the radiative phase sets in, the bubble pressure drops and no longer drives the shock 
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into the ISM. The expansion rate of the shell then slows significantly following an expansion 
law of r ~ t 1 / 4 until the subsequent wind momentum exceeds the present momentum in the 
shell. Radiative losses are important when t = Tcool, defined as the time when the total 
energy radiated by the bubble is equal to the total energy input from the wind. The cooling 



time Tcool is given by the following equation based on that given by iMcCray fc Kafatos 



(119871 ) (for solar metallicity), 



COOL 



7.3 x 10- 5 L^ 3 n- a7 yr (5) 



An additional effect of bubble cooling is that the acceleration of relativistic particles will 
shut off at t = Tcool, since the bubble shock collapses to the outer radiative shell where the 
higher gas densities preclude efficient particle acceleration. The low acceleration efficiency is a 
consequence of excessive collisional losses at the low energy phase of the acceleration process. 
Accordingly, this test requires that the bubble be in the adiabatic phase, i.e. Tcool > t. 



3.4. Pressure confinement 



The expanding shell is subject to confinement due to the outside ambient gas and mag- 
netic pressures of the ISM. This pressure stops the expansion w hen it becomes comparable 
to the ram pressure of the expanding shell (e.g. lOey et al.ll2002l ). We define the confinement 
time Tconf to be the time at which the shell expansion slows to lOkms -1 , adopted as the 
effective sound speed of the ISM in M82. We thus require Tconf > t, where the expansion 
rate is given by dr/dt = I r/t. 



3.5. Presence of a detectable HII region surrounding the cluster 

The most massive stars in the cluster will produce a surrounding HII region which 
should be detectable as a thermal radio source. Such emission is probably the source of the 
compact HII regions seen in M82, which so far have no detectable non-thermal emission. 
Since the expanding shell should be either ionization bounded, or nearly so (see §5.2.3), 
the expected thermal or free-free flux density Sexp is derivable directly from the Lyman 
continuum luminosity Ly given by equation (3e). The expected (optically thin) radio flux is, 

S EX P = 1-32 z.- 01 Ly 5 oD- 2 mJy, (6) 

where Lyso is in units of 10 50 photons s _1 , and D is the distance in Mpc (in this case 3.25 Mpc). 
This emission is not significantly attenuated by dust, and is characterized by its flat or rising 
radio spectrum, compared to the non-thermal emission which has a negative slope. 
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The difficulty with applying this test to the SNR sources is the confusion by non-thermal 
emission. Since the spectrum of optically thin free-free emission is flat, detectable thermal 
emission would manifest itself as a positive curvature of the radio spectrum, affecting the 
spectrum preferentially at high frequency (in this case 15-20 GHz). In order to search for 
a thermal component, the radio spectrum for each member of the input sample was fit by 
a curve representing combined thermal and non-thermal emission. In no case was there 
evidence for a thermal component. A 3<j upper limit Sobs was therefore assigned in each 
case, where a is the standard error of the thermal component derived from the fit. Thus 
the test in this case requires that the expected flux not exceed the observed upper limit, i.e. 
Sexp < Sobs j where the typical value for Sobs is a f ew tenths of a mJy 



4. RESULTS 

In applying each of the limits above, a violation by a factor of two was permitted to 
allow for the uncertainty in deriving the physical parameters. For example, the energy test 
in §3.1 was regarded as successful as long as the energy requirement of relativistic electrons 
did not require the parameter f to exceed the prescribed value of 0.1 by more than a factor 
of two. 

For each source in the input sample, the five tests were applied considering ranges 
in possible SSC mass and ISM density of 10 3 - 1O 6 M and 10 - 10 5 cm~ 3 , respectively. 
The logarithmic step size was 0.5 for both parameters. The result is an output sample 
comprising sources which passed all five tests described in §3 for a specified range in (M4, 
no). Briefly stated, the mass range is constrained by tests 3.1 and 3.5, and for a given 
mass (i.e. wind power), the ISM density is constrained on the low side by ages too low to 
provide sufficient energy in test 3.1, and by tests 3.2 and 3.3 on the high side by excessive 
radiation or cooling losses. Test 3.4 did not provide any limiting constraint within the ranges 
allowed by the other tests, signifying that the shells are not pressure confined. Table 1 lists 
the members of the output sample with the inferred medians and ranges of SSC mass and 
ISM density. The output sample contains all but two members of the input sample, namely 
42.5 + 61.9 and 42.8 + 61.3, both of which failed test 3.5. Failure of this test signifies that the 
masses (and hence Lyman continuum luminosities) required were so high that the predicted 
free-free emission from the associated HII regions was higher than the upper limits set by 
observations. Thus it may be concluded that the WDB model is consistent with the great 
majority of CNRs, though this does not preclude their interpretation as SNRs. It remains 
distinctly possible as well that CNRs comprise a mixture of both types. 

Figure 2 shows the distributions of the median values in Table 1. The distribution shows 
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no significant correlation between M4 and no, and the centroid of this distribution corresponds 
to M = 1.8 x 10 4 M^ and n n = 2.6 x 10 3 cm -3 . Since the mass range for observed SSCs 



is about lO 5 ^ 1 ) flMelo et al.1 120051 ). we conclude that if the CNRs are WDBs, the SSCs 
responsible are at the low end of the mass range for M 82. The derived value of n is consistent 
with the density associated with the inter-clump medi um of M82 (~ 10 3 c m" 3 ) rather tha n 
the molecular cloud densities (~ 10 5 cm -3 ) (see for e.g. iGiisten et al.lll993l ; iMao et al.ll2000l ). 
which suggests that the clusters would have consumed most of the dense molecular gas in 
their parent clouds. 

It must be emphasized that the tests do not lead to a discrimination between the 
WDB and SNR models. In particular, some sources in the output sample exhibit variability 
consistent with the SNR hypothesis, and a small number have measured expansion velocities 
of several thousand kms" 1 suggesting an SNR in free expansion. However, it is interesting 
that at least two of sources in the latter category show no measurable decre ase in flux as 



would be expected for a freely expanding cloud of particles and magnetic field. iFenech et al. 
( 120051 ) report an expansion velocity of 10,500 kms -1 for 43.2 + 58.4, based on an apparent 



increase in size of 0.7 mas per year between 1992 and 2002. This figure corresponds to an 
increase of about 1 % per year in radius. H owever, a linear least squares fit to the data for 
radio flux vs. data in lKronberg et al. yields a chang e in flux of +0.50±0.54% per year 

over the period 1982-1992. Similarly. iBeswick et al.l (120061 ) derive an expansion velocity near 
10,000 kms" 1 for the source 43.3 + 59.2. Their Figure 8 shows that the rate of increase in 
source size, extrapolated back to 1982-1992 yields at least 4% per year, whereas the flux data 
yield a change of +0.40 ± 0.22% per year over this period. These apparent inconsistencies 
suggest that it may yet be premature to use the measured expansion velocities to conclude 
anything definitive about the nature of these sources. 



5. DISCUSSION 

5.1. Possible similar origin between CNRs and compact HII sources 

Though there is no known relationship between the non-thermal and thermal compact 
radio sources (i.e. compact HII regions), it is plausible that both m ay be associated with 



SSCs. The median radius of the compact HII regions measured by iRodriguez-Rico et al. 



( 120041 ) is about 4pc, comparable to and somewhat larger than those for the CNR population 
(lpc). The HII regions may be excited by massive stars in the clusters, which should also 
produce WDBs. Proceeding on the assumption that the comp act HII regions are associate d 



with SSCs, we used their 43 GHz flux densities measured by IRodriguez-Rico et al.l (120041 ) 



together with equations (3e) and (6), to derive the masses of the associated clusters, assuming 
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Fig. 2. — Plot in 3-D showing the range in parameter space for cluster mass and ambient 
ISM density possible for a WDB interpretation of a representative sample of compact non- 
thermal radio sources. The individual values are medians of the range of acceptable values 
for each source (see Table 1). 

that they have the same IMF as those associated with the CNRs. Figures 3(a) and 3(b) 
show a comparison between the inferred distribution of cluster masses underlying CNRs 
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(from Figure 2), and the corresponding distribution for the HII sources. The compact HII 
regions are seen to be associated with significantly more massive SSCs. The median mass for 
the HII region clusters is 1.1 x 1O 5 M compared to 1.8 x 1O 4 M for the CNR clusters. By 
the WDB hypothesis then, the compact HII regions would be associated with young clusters 
of mass typical for M82, while the CNR clusters would be under-massive. 

If compact HII regions are associated with cluster WDBs, and if WDBs invariably 
produce synchrotron emitting electrons, then the lack of detectable non-thermal emission 
from the HII sources is puzzling. The absence is possibly attributable to shorter inverse 
Compton lifetimes associated with the higher cluster luminosities. These bubbles may also 
be older, in their radiative phase, and/or pressure confined. In this case particle acceleration 
will have ceased. 

In order to identify the most likely possibilites, we adopt M = 1.1 x 1O 5 M for the 
underlying clusters, r = 4 pc for the associated bubble radius, as described earlier, and 
H = 10~ 3 gauss for the emitting region. We applied the same tests to this case as for the 
CNRs, except for test 3.5, since the free-free flux in this case was used to estimate the cluster 
mass. We find that nonthermal emission is possible provided n < 10 4 cm~ 3 , for which the 
bubble age inferred from equations (1) and (3b) is less than about 10 5 yr. At higher den- 
sity the corresponding age would be greater than the inverse Compton lifetime for electrons 
emitting at 5 GHz. At no ~ 10 5 cm -3 radiative cooling also becomes an important factor, 
shutting off the supply of relativistic particles. The conclusion is that synchrotron emission 
would be supressed in the observed HII sources if the ambient gas density exceeds 10 4 cm~ 3 . 



5.2. Properties of the WDB shells for the CNRs 

5.2.1. Mass, column density and age 

Figures 4, 5, and 6 show the distributions of the inferred mass, column density and 
age for the shells in the WDB hypothesis. These parameters were derived from the median 
values of (M, n ) for each member of the output sample in Table 1. The shell masses are 
derived from the source volume and the ISM density. The column density of hydrogen N s 
associated with the shell thickness is derived using Ns = 1/3 n r, appropriate for a thin shell, 
and the ages are derived from equation (1). 

Note that the shapes of the distributions reflect largely the uncertainties in the determi- 
nation of these parameters rather than their intrinsic distributions. It is the median values 
that are the most relevant. The median shell mass and column density are 212 M and 
2.3 x 10 21 cm~ 2 , respectively. Using a figure of 1.87 x 10 21 cm -2 for a visual extinction of 
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Table 1: Physical parameters, derived SSC masses, and interstellar medium densities for an 
acceptable fit of the CNRs to the WDB model for the output sample 



Name 


d a 


Q b 






Hi c 




logW 




log n g 




pc 


mJy 




10 48 erg 


10~ 4 gauss 




(M 


) 


( 


cm -3 ) 














(median, 


range) 


(median, range) 


39.1+57.4 


1.7 


8.6 


-0.38 


3.3 


7.6 


4.3 


(4.0, 


4.5) 


3.4 


(2.5, 4.0) 


39.4+56.1 


3.3 


3.1 


-0.21 


5.0 


3.4 


4.2 


(3.5, 


4.5) 


2.9 


(1.5, 4.0) 


39.6+53.4 


2.8 


4.4 


-0.71 


4.8 


4.3 


4.0 


(4.0, 


4.0) 


3.5 


(3.0, 4.0) 


39.8+56.9 


3.1 h 


2.1 


-0.49 


3.4 


3.0 


4.2 


(3.5, 


4.5) 


2.7 


(1.0, 4.0) 


40.3+55.1 


1.6 


2.3 


-0.55 


1.4 


5.5 


4.2 


(3.5, 


4.5) 


3.1 


(1.5, 4.0) 


40.6+56.1 


2.9 


3.8 


-0.72 


4.7 


3.9 


3.9 


(3.5, 


4.0) 


3.3 


(2.5, 4.0) 


40.7+55.1 


2.0 


7.9 


-0.52 


5.9 


8.0 


4.5 


(4.4, 


4.5) 


3.2 


(3.0, 3.5) 


41.3+59.6 


1.1 


8.6 


-0.54 


2.0 


11.1 


4.3 


(4.4, 


4.5) 


3.4 


(2.5, 4.0) 


42.7+55.7 


4.1 


4.4 


-0.61 


7.1 


2.9 


3.5 


(3.5, 


3.5) 


3.5 


(3.5, 3.5) 


43.2+58.4 


1.0 


15.3 


-0.67 


2.6 


14.5 


4.5 


(4.5, 


4.5) 


3.7 


(3.5, 4.0) 


43.3+59.2 


0.6 


30.3 


-0.64 


1.9 


28.0 


4.5 


(4.5, 


4.5) 


4.0 


(4.0, 4.0) 


44.0+59.6 


0.8 


62.0 


-0.51 


3.8 


25.6 


5.0 


(5.0, 


5.0) 


3.5 


(3.5, 3.5) 


44.3+59.3 


1.9 


6.7 


-0.56 


3.3 


6.3 


4.3 


(4.0, 


4.5) 


3.2 


(2.0, 4.0) 


44.5+58.2 


2.2 


7.2 


-0.61 


4.2 


5.8 


4.0 


(4.0, 


4.0) 


3.5 


(3.0, 4.0) 


44.9+61.1 


1.6 


6.6 


-0.45 


2.6 


7.3 


4.0 


(4.0, 


4.0) 


3.7 


(3.5, 4.0) 


45.2+61.2 


1.1 


24.1 


-0.68 


4.1 


15.4 


4.5 


(4.5, 


4.5) 


3.5 


(3.5, 3.5) 


45.3+65.2 


2.0 


5.4 


-0.62 


3.1 


5.7 


3.9 


(3.5, 


4.0) 


3.6 


(3.0, 4.0) 


45.4+67.4 


2.2 


4.7 


-0.57 


3.5 


5.2 


4.0 


(4.0, 


4.0) 


3.5 


(3.0, 4.0) 


45.5+64.8 


3.1 h 


1.7 


-0.15 


3.8 


3.2 


4.2 


(3.5, 


4.5) 


3.0 


(1.5, 4.0) 


45.8+65.3 


2.1 


5.1 


-0.23 


3.8 


5.8 


4.3 


(4.0, 


4.5) 


3.2 


(2.0, 4.0) 


45.9+63.9 


2.2 


4.1 


-0.53 


3.2 


4.9 


4.3 


(4.0, 


4.5) 


3.2 


(2.0, 4.0) 


46.5+63.9 


1.5 


5.2 


-0.57 


2.0 


7.3 


3.9 


(3.5, 


4.0) 


3.6 


(3.0, 4.0) 


46.6+73.8 


3.1 h 


3.7 


-0.78 


5.8 


3.9 


4.0 


(4.0, 


4.0) 


3.2 


(3.0, 3.5) 


46.7+67.0 


3.0 


5.2 


-0.57 


5.0 


4.0 


4.3 


(4.0, 


4.5) 


3.2 


(2.0, 4.0) 


47.4+68.0 


3.1 h 


4.1 


-0.83 


6.6 


4.2 


4.3 


(4.0, 


4.5) 


2.9 


(2.0, 3.5) 



"source diameter 

'flux at 1 GHz 

c spectral index (S = So v a ) 

d equipartition energy = Ei x (1 + k) 4 / 7 </> 3 / 7 

e equipartition magnetic field = Hi x [(1 + k)/(/>] 2 / 7 

^SSC mass (median and range) 

9 ISM hydrogen density (median and range) 

^based on assumed diameter of 0.2" 



Note. — Source d ata obtained from iHuang et all ( 19941 ). iMuxlow et all ( 19941 ). lAllenl ( 1999h and 
McDonald et all |2002l ) 
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Fig. 3. — Distribution of (a) allowed masses for the WDB model of the compact non-thermal 
sources, based on Figure 2 and (b) the same for comp act HII regions based on the observed 
43 GHz flux densities by iRodriguez-Rico et al.l (120041 ) . 
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1 magnitude (jBohlin. Savage fe Drake! Il978l ). the median shell extinction is 1.2 m , indicating 
that many of the shells in this picture are partially transparent to optical radiation. Thus 
any neutral gas would be significantly photo-dissociated, comprising mainly HI instead of 
H 2 . The median age of the winds is about 1.3 x 10 4 yr, significantly greater than for the 
SNR hypothesis 10 3 yr). The derived ages have implications for the star formation rate 
in clusters with masses near 10 4 M Q . From the median mass of ~ 2 x 1O 4 M for the 25 
SSCs in the output sample and adopting a maximum age of ~ 5 x 10 4 yr, the implied SFR 
is ~ 10M Q yr _1 . This value seems not unreasonable compared to recent estimates for the 
total SFR of up to lO-SSMoyr" 1 jF5rster Schreiber et alJliooi l. 




Fig. 4. — The distribution of log Mass contained in the shells swept up by the expanding 
bubble in the WDB interpretation of compact non-thermal sources, computed from the size 
and median gas density associated with each source in Table 1. 
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log Column density (cm" ) 



Fig. 5. — The distribution of log Column Density associated with the shells swept up by the 
expanding bubble in the WDB interpretation of compact non-thermal sources, computed 
from the size and median gas density associated with each source in Table 1. 

5.2.2. Magnetic field and gas density 

The magnetic field is not directly observable, but is assumed to be the equipartition field 
derived from the synchrotron luminosity. The median value is 5.7 x 10~ 4 (1 + k) 2 / 7 $~ 2 / 7 
gauss. However, the WDB model provides an independent (though model dependent) 
method for estimating the magnetic field in the expanding shell using pressure balance 
between the field in the shell and the ram pressure of the ISM. Assuming for the moment 
that gas pressure in the shell is not a significant contributor, and that the magnetic field is 
the only means of support, the field strength is given by H 2 /8 n p Vg where V s is the shell 
expansion velocity. Taking Vs = | r/t, the median shell speed for the sources is ~56kms _1 , 
and the corresponding value for the magnetic field in the shell is 1.5 x 10 -3 gauss. This 
figure is close to median equipartition field if allowance is made for the factor (l+k) 2 / 7 $~ 2 / 7 
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Fig. 6. — The distribution of log Age associated with the shells swept up by the expanding 
bubble in the WDB interpretation of compact non-thermal sources, computed from the 
median mass (and associated wind power) and median gas density associated with each 
source in Table 1. 



which will exceed unity. It is interesting to compare this value with t he ambient inter- 
stellar field in M82, which is b e tween 5 x 10~ 5 and 1.4 x 1CT 4 gaus s (IRieke et al.1 Il980l ; 



Klein. Wielebinski fc Morsil Il988l ; ICondonl 1 1992 ; iThompson et al.l 120061 ) . Assuming a value 
of 10~ 4 gauss for the interstellar field and adopting 1.5 x 10~ 3 gauss for the field in the 
shell, we find that interstellar field is compressed by about one order of magnitude. Assum- 
ing the compression is 1-D with the field parallel to the shock, and that the field is frozen 
in, the shell gas density n s corresponding to the median value of no = 2.6 x 10 3 cm~ 3 is 
n s = 3.9 x 10 4 cm~ 3 . Two factors which could act to lower this density are (1) a lower value 
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for the shell magnetic field by taking account of other possible sources of pressure, and (2) a 
higher value for the i nterstellar magnet i c field , i.e. stronger than the equipartition value, as 
suggested recently by lThompson et all (120061 ). Since this density cannot provide significant 
pressure at the low temperatures expected for the radiative shell, the initial assumption of 
ignoring gas pressure appears justified. 



5.2.3. Fractional ionization 

The hydrogen emission measure (EMb) of a geometrically thin ionization bounded shell 
can be determined directly from the cluster Lyman continuum luminosity, which when com- 
bined with equation (3e), yields 

EM B = 4.3 x 10 6 M 4 /r 2 cm" 6 pc (7) 

This equation assumes an ionized gas kinetic temperature of 10 4 K and Menzel case B recom- 
bination. For the median cluster mass of 1.8 x 10 4 M , equation (7) yields EMb = 7.7 x 10 6 
cm~ 6 pc, using a mean radius of 1.0 pc. This value may be compared with that occurring if 
the shell were fully ionized, which is, 

EM S = 3.2 x 10" 19 n s N s cm" 6 pc, (8) 

where Ng is the column density in cm -2 . With n s = 3.9 x 10 4 cm -3 and Ns = 2.3 x 10 21 cm" 2 
(see § 5.2.1 and 5.2.2), equation (8) yields EM S = 2.9 x 10 7 cm" 6 pc. Thus the ratio EM B /EM S ~ 
0.3, indicating that a typical shell is probably ionization bounded. However, the ratio is so 
close to unity that some shells would be substantially neutral and others substantially ion- 
ized. 



6. ENERGETICS FOR WDB AND SNR MODELS 

The radio luminosities of the CNRs are higher than for galactic SNRs. The total kinetic 
energy for SNR ejecta is 10 51 erg, whereas the total energy available in the WDB model for 
the median mass and source age is about 4 x 10 49 erg. Thus at first sight, the efficiency 
of relativistic particle acceleration in the WDB interpretation needs to be at least an order 
of magnitude higher than for the SNR interpretation. We examine this argument in more 
detail below. 
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6.1. Relativistic particle efRciencies in Galactic SNRs 



Some estimates of the efficiency of relativistic particle acceleration rj in SNRs are 
based on a comparison between the total equipartition energy for the relativistic parti- 
cles and__the__total_jdne_Uce of 10 51 erg for the ejecta. An early estimate of this type 
by iGinzburg & Syrovatskiil (119641 ) led to rj — 1 — 10%, but with very large uncertainties. 
A primary difficulty with this method lies in the uncertain ratio k of energy between pro- 
tons (which radiate inefficiently) and electrons (which provide essentially all of the observed 
synchrotron emission). Another method is based on direct estimates of the local cosmic 
ray energy density and the assumption that these are generated solely by SNRs. For this 
method, estimates of the cosmic ray residence time in the disk and the rate of Galactic 
SNe are required. These estimat es vary widely becau s e of t he uncertainty in the parameters 
used. Using the la t ter m ethod, iBlandford &: Eichlerl (119871 ) find rj = 3%. At the other ex- 
treme, iFields et al. (booih . using the same method, find rj = 30%, admitting an uncertainty 
of a factor of at least three. Theoretical e stimates also tend to be rather high, with fig- 



ures that reach 50% (IBerezhko Volkl 119971 ) . Note that in local cosmic rays the fraction of 
electrons is only ~ 2%, leading to corresponding electron acceleration efficiencies r] e ranging 
from -0.06% to -1%. 

Recent studies of TeV gamma ray emission produced by inverse Compton scattering 
of cosmic background photons by energetic electrons provide fresh insight into the effi- 
ciency of electron a cceleration, concluding that in some SNRs rj e is in the range 1.0-2.5% 
(IKeshet et al.ll2003l . and references therein). The corresponding fraction of shock thermal 
energy carried by relativistic electrons is then about 5%. Note that the proton/electron 
energy ratio k must then be much lower than 50, and p erhaps nearer 10 to keep the proton 
contribution from dominating the thermal gas pressure. Ellison et al.l (120041 ) estimates that 
the fraction of post shock thermal energy occupied by relativistic particles (including ions) 
in SNR shocks is about 50%. If we regard this as a reasonable upper limit, then combining 
this with the figure of 5% for the partial pressure of electrons, yields k ^ 10. In order to be 
compatible with the observed radio luminosities of Galactic SNRs, the large fraction of shock 
thermal energy in the form of relativistic particles implies a high ratio of particle /magne tic 



field energy, so that at least some SNRs are far from equipartition (e.g. iDyer et al.l 1200 ll ). 



Presuming that the non-thermal sources in M82 have a WDB origin with relativistic 
particle efficiency comparable to values recently determined for galactic SNRs, but with 
higher (equipartition) magnetic fields, then a basis exists for the high ratio of radio luminosity 
to mechanical input energy for the WDB sources compared to galactic SNRs. Compared to 
Galactic SNRs, the higher ratio of radio luminosity to available particle energy in WDBs 
in M 82 is attributable to its strong interstellar field (1-2 orders stronger than the Galactic 



field), further amplified by the shocks associated with the expanding shells of the WDBs. 



6.2. Relativistic particle efficiency in the WDB model 

In §3.1, the ratio f = Ei/E w was set to f = 0.1 for the conservation of energy test. Here 
we further justify this choice by showing that the associated electron acceleration efficiency 
is comparable to that estimated for Galactic SNRs. 

The efficiency of converting wind mechanical energy to relativistic particle energy is 
given by: 

E p /E w = if(l + k) 4 / 7 $ 3 / 7 , (9) 

where the factor | accounts for equal amounts of Ei going into particle and magnetic energy. 
Since the electron energy is given by E p /(l+k), the efficiency of conversion to relativistic 
electron energy rj e is 

Ve = \f[<S>/(l+k)f 7 (10) 

Setting f = 0.1 and adopting 1 < (1 + k) < 10, we obtain 0.02$ 3 / 7 ^ r? ^ 0.05 $ 3/7 , com- 
parable to that in at least some Galactic SNRs. Thus f = 0.1 is seen to be a reasonable 
choice based on comparisons with recent results for SNRs. 



6.3. Expansion, ionization and bremsstrahlung losses 

No account was taken of expansion, ionization or bremsstrahlung losses by the rela- 
tivistic particles in the tests described in §3, chiefly because they are less certain and/or 
less stringent. Expansion losses may be incurred as the particles propagate outward from 
the inner bubble shock toward the expanding outer shock, where the observed radiation is 
assumed to occur. Adiabatic losses scale inversely with (V2/V1) 1 / 3 where Vi and V2 are the 
initial and final volumes occupied by the particles. We examine the worst case in which Vi 
and V 2 constitute the entire volumes bounded by the inner wind and outer shock radii R sw 
and R s , re spectively. These rad ii vary as R sw ~ t^/^ and R s ~ t 3// ^, as given by equations 4.2 



and 4.4 of iKoo fc McKed (119921 ) . For conditions relevant to this paper, these equations yield 
a current value of (R sw /R s ) ~ 1/4. Using these expansion laws, assuming a steady injection 
of particles into the bubble at fixed fractional radius R sw , and integrating over R sw since 
t=0, we estimate that the total energy loss is about 80% for a cluster with a current value 
of R sw /R s ~ 1/4. However, we emphasize here that this is a pessimistic estimate, and also 
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a quite uncertain one especially since the geometrical factors adopted here would be quite 
different when the source of the wind has a finite radius. 

Relativistic electrons co-existing with gas in t he de nse shell will suffer collisional (i.e. 



ionization) losses. Using equation 4.411 in iLangl (119981 ) for the ionization loss, it may be 



shown that electrons emitting at frequency z/5 (z//5GHz) in a medium with neutral density 
n4 (n/10 4 cm~ 3 ) and magnetic field H_3 (H/10 -3 gauss) lose energy on a time scale 

Ti = 1.0 x 10 4 nj 1 HlJ /2 z/ 5 1/2 yr. (11) 

Adopting n 4 = n s = 3.9, H_ 3 = 1.5 for the shell, as suggested in §5.2.2, then Ti ~ 2 x 10 3 
years for electrons radiating at 5 GHz, which is less than the median age of the bubbles. 
Because of the energy (or frequency) dependence of this timescale, the effects of such losses 
would be most severe below 5 GHz, and would produce a spectral turnover or cutoff at low 
frequencies. In fact, a turnover occurs below 1 GHz in most sources in M 82, but is attributed 
to the effects of free-free absorption by ionized gas in the ISM. 

A similar argument for the collision loss can be made if the shell gas is ionized, and the 
collision loss timescale (computation not shown here) in this case is estimated to be about 
a factor two shorter than the value for neutral gas. 

If the shell gas is ionized, then t he rela t ivistic electrons will also lose energy by bremsstrahlung 



radiation. From equation (4.417) of ILangl (119981 ). it may be shown that the loss time is only 



logarithmically dependent on particle energy (or frequency), and is given approximately by 

T B = 5.4 x lf^njV (12) 
for electrons radiating near 5 GHz in a field of 1.5 x 10~ 3 gauss. 

Again adopting n 4 = n s = 3.9, Tb = 1.4 x 10 3 yr, once again less than the median 
age of the bubbles. Thus it appears that ionization and/or bremsstrahlung losses would be 
significant if the relativistic particles were confined to the dense shell. However, it seems 
unlikely that the relativistic particles would be confined strictly to the dense shell, but would 
occupy a somewhat larger volume, likely encompassing the thin shell. 



THE SIZES OF CNRs AND SSCs 



The diameters of the CNRs in the output sample range from ab out 0.6 to 4. 1 pc. T he 



mean photometric diameter for SSCs in M 82 clusters observed by iMelo et al.l (120051 ) is 
d = 11.4 ± 2.8 pc, largely independent of cluster mass. This difference clearly raises ques- 
tions about the WDB picture since the clusters ought to be smaller than the bubbles they 
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produce. However there are two considerations which could in principle resolve this difficulty. 
One is a possibility that mass segregation has occurred, with the most massive stars now at 
the center of the cluster. Although some observational evidence exists for this from color 
gradients (IMcCradv. Graham fc Vaccall2005f) . the matter is unresolved, since there is also con- 



tradictory evidence ( Smith et al. 20061 ) . Secondly, the hypothesized clusters underlying the 
CNRs are extremely young, since their winds are ~ 10 4 years old. Such clusters may be more 
compact than older clusters, which are likely modified by tidal interactions. Applying the 
Virial Theorem to a 1O 4 M cluster with a velocity dispersion of lOkms -1 , comparable to gas 
turbulent velocities within molecular clouds, yields d~0.5pc, smaller than those observed 
by more than an order of magnitude, and smaller than the CNR diameters. Thus a possible 
explanation is that the youngest clusters have large binding energies, which decrease with 
age due to energy input from tidal interactions with molecular clouds and other clusters. 



8. THE CNRs IN OTHER WAVEBANDS 

M82 has been mapped in a number of wavebands from the radio to X-Rays, and such 
images could in principle yield evidence bearing on the relevance of the WDB origin for the 
CNRs. 



8.1. Near-IR 



Alonso-Herrero et al.l (120031 ) compared HST images in the [Fell] 1.644 /im line and ad- 
jacent continuum with the locations of 44 CNRs at sub arcsec resolution. The [Fell] lines are 
collisionally excited, and are well known signatures of SNRs both in our Galaxy and nearby 
galaxies. They find that 30 - 50% of the CNRs possess associated line emission, which they 
consider support for the SNR interpretation. The observed luminosities (corrected for ex- 
tinction) are in the range 3 x 10 36 — 2.2 x 10 38 ergs _1 . 

For comparison, we have estimated the [Fell] luminosity expected for a typical WDB 
shell. There are two possible sources - one with luminosity L s , associated with the outer 
shock, and the other with luminosity L p , associated with gas in the shell photoionized by 
the cluster. It is shown in the Appendix that these may be written as 

L s = 2.0 x 10 37 ^ S M 4 ergs" 1 (13) 
L p = 3.3 x 10 37 77 P M 4 erg s" 1 (14) 
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where rj s and 7] p are the ratios of [ Fell] l uminosity to shock kinetic energy and Paschen/? 
luminosity, respectively. iMouri et al.l (120001 ) estimate rj s = 2 x 10~ 3 for M82. For rj p , their 
Figure 8(a) shows a plot of rj p vs. ionization parameter U for a limited set of conditions. 
In our case U = Ly/4 7rr 2 n s c, where c is the speed of light. Assuming solar abundances, 
n s = 10 4 cm -3 , and ionization by a blackbody with T=40,000K, then for M4 = 2.0, and 
r = 1 pc, we estimate f] p = 0.3. The resulting luminosities are L s = 8.0 x 10 34 ergs _1 and 
L p = 2.0 x 10 37 ergs _1 , showing that the emission from the photoionized shell dominates 
that from the shock by more than two orders of magnitude. Comparison with the observed 
luminosity range above shows that the [Fell] luminosity of a typical WDB would be in the 
middle of the observed range. 

The ratio rj p may also be used to predict the ratio [F eII]/Br 7 for comparison with the 
observed ratios in Figure 10 of Alonso-Herrero et al. J2003h . The predicted ratio is about 0.4, 
compared to a median observed value near 1.0. Both are less than that shown for Galactic 
SNRs which is about 10. The predicted values of L p and [Fell] /Br 7 would be higher for lower 
ionization parameter, for e.g. if the shell density were as high as our estimate of 4 x 10 4 cm -3 
discussed in section 5.2.2. Higher densities still would not help however since this value is 
approximately the critical density for this transition, setting a corresponding limit on L p . 



It is also shown in IMouri et al.l (120001 ) that if the abundances are sub-solar, with 90% 
of Fe locked up in grains, then rj p and hence L p would be lower than the above estimate by 
an order of magnitude, and a t ypical WDB cluster probably would not have been detected 
by I Alonso-Herrero et al.l (120031 ) . 



Alonso-Herrero et al.l (120031 ) did not conduct a search for coincidences between CNRs 



and continuum near-IR sources. Their Figure 2(a) shows the image in the adjacent con- 
tinuum band at 1.64 /im. Using a position and flux calibrated electronic version of this 
image provided by the authors, we searched for coincidences between 1.64 /im continuum 
point sources and both CNRs and thermal compact sources. For the CNRs two coincidences 
were found with separation ^0.5", which is consistent with chance associations (estimated 
probability 20-30%). For the thermal sources there were two coincidences with separation 
^0.2", which may also be chance associations. Thus there is no significant indication that 
the near-IR continuum and radio populations are physically associated. 

We estimate that the faintest detectable continuum point source on these images is 
about 0.3 mJy. For comparison, SB99 predicts that, with cluster parameters used in this 
paper, an expected 1.64//m continuum flux of 0.12M 4 mJy from the cluster stars. With 
an estimate of typically 2 m of extinction at H (based on a typical value of 1.4 m at K-band 
estimated by Alonso-Herrero et al.), we conclude that if there are clusters associated with 
compact radio sources, they would likely be below the detection threshold of these images. 
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8.2. Mid-IR 

IR emission at wavelengths near 10 /im is expected in the WDB model from dust in the 
shells heated by optical radiation from the underlying clusters. The shells swept up by the 
shock are nearly opaque in the visual region as shown in §5.2.1, so it may be assumed for 
present purposes that all of the bolometric luminosity is converte d to IR emission. Recent 



high resolution (0.4") observations of the central 400 pc of M82 by iLipscy fc Plavchanl (120041 ) 
at 11.7 and 17.65 fim reveal seven star-forming clusters which together provide ~ 15% of the 
total mid-IR luminosity. The authors kindly provided electronic versions of these images to 
permit a search for associations with the CNRs. A search at both wavelengths was made and 
failed to yield any coincidences with the output sample members. In general, the upper limits 
for coincident sources were about 50mJy and 100 mJy at 11.7 and 17.65 /zm respectively. 

A crude and optimistic estimate of the expected level of mid-IR flux was made by 
assuming, as in the case of the Milky Way, that small grains absorb 20-30% of the incident 
optical light a nd re-radiate this energy in the mid-IR as a result of impulsive grain heating 



(iDraind 120031 ) . The remaining optical light is radiated in the far-IR by larger grains under 
equilibrium conditions. By assuming that the mid-IR emission acts like a black body with 
peak emission at 12 /im and total luminosity equivalent to 30% of the bolometric luminosity 
given by equation (3d), we obtain 44 mJy for the flux associated with the 2 x 10 4 M cluster, 
just below the mean upper limit of 50mJy at 11.7 //m. This result suggests that such clusters 
might be detectable if the sensitivity were improved by a factor of several. 



8.3. X-Ray emission 



Emission at soft X-Rays is expected from the hot expanding bubble. High angular 
resolution observations of M82 using Chandr a show about 25 discrete sources down to a 
luminosity detection limit of 2 x 10 37 ergs _1 (IZezas et al.l 120041 ) . Apart from the powerful 
variable source tentatively identified with the CNR 41.30 + 59.6, these sources appear in the 
majority to be X-Ray binaries associated with young star clusters, and not identified with 
CNRs. 



To estimate the expected X-Ray luminosity of a typical CNR in the WDB picture, we 
first note that the total wind power for a 2 x 1O 4 M cluster is 1.1 x 10 38 ergs _1 which 
represents an upper limit to the X-ray luminosity of the bubble. The true X-Ray luminosity 
would be only a fraction of this amount beca use: (1) the thermal energ y of the bubble is 
5/11 that of the total wind mechanical energy (IMac Low feMcCrayl 119881 ); (2) the bubble is 
in the adiabatic expansion phase and the cooling rate at this stage is lower than the wind 
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power; and (3) most of the cooling radiation is emitted in the interface regio n between the 



bubb le shock and outer shell by UV fluorescence lines such as OVI A 1035 (IWeaver et al. 



19771 ). Consequently, the X-ray luminosity of a typical bubble could be weaker than the 



current detection limit. 



8.4. Association with molecular clouds 



Using CO J=2-l interferometric observations of M 82 at sub-arcsec resolution, iKeto. HofcLo 



(120051 ) resolved the molecular gas in the starburst region into a large number of compact 
clouds with masses ranging from 2 x 10 3 to 2 x 10 6 M Q . Surprisingly, the CNRs seem 
to cluster around the periphery of molecular emission regions, and there is almost an anti- 
correlation between molecular clouds and the locations of the compact radio sources. The 
authors suggest that, since the CNRs must be signposts of recent star formation, the as- 
sociated gas has either been fully consumed by stars, or any associated gas has a spatial 
scale > 10". Perhaps it is noteworthy that in §4 it was shown that in the WDB model, 
the ISM into which the bubbles are expanding is about 2.6 x 10 3 cm~ 3 , comparable to that 
estimated for the warm diffuse inter-clump medium (~ 10 3 cm~ 3 ) rather than the cooler and 
denser clumps (~ 10 5 cm -3 ) where stars form. Thus the clusters could have consumed the 
dense gas, and are now expanding into a remnant warm, primarily dissociated medium, 
invisible in CO emission. 



9. SUMMARY AND CONCLUSIONS 

We have introduced and tested the hypothesis that many of the CNRs in M82 are 
associated with WDBs rather than SNRs. A primary motivation for this hypothesis is 
that most of the radio sources exhibit no discernible time variation in their flux densities, 
suggesting an age greater than expected for SNRs, but comfortably within the range expected 
for WDBs. In the WDB picture, the synchrotron emitting particles are accelerated by the 
first order Fermi mechanism at the shock at the site of the interaction between the cluster 
wind and bubble gas. 

The WDB model was tested using a limited but representative sample of objects with 
well defined radio spectra. The model, based on the theory of WDBs and model parameters 
from SB99, was constrained to satisfy the energy requirements of the synchrotron sources 
and limits imposed by losses of the synchrotron emitting particles and radiative cooling time 
scale of the bubble. Further constraints were that the predicted free-free flux density from 
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the HII region produced by the young cluster not exceed the upper limits set by observation, 
and that the bubble age not exceed the time to pressure confinement. These considerations 
constrain the associated cluster masses to ~ 2 x lO 4 ^ 1 - ) M Q and the surrounding ISM density 
to ~ 3 x lO 3 ^ 1 ) cm -3 . The corresponding ages cover the range from a few x 10 3 years to a 
few x 10 4 yr. This range in age is shown to be consistent with a plausible star formation rate 
in M 82. Estimates of the flux densities in various wavebands of clusters at the median mass 
of 2 x 10 4 M show that published observations are not yet sensitive enough to confirm (or 
rule out) such clusters, but in some cases are very close. 

Thus, we conclude that the WDB hypothesis appears to be a plausible alternative to 
the SNR picture for the origin of the CNRs, but the tests do not rule out an SNR origin. 
It is also possible that the CNRs comprise a mixture of both types. It is clear that some 
critical tests are needed to discriminate between the two ideas. The most stringent tests 
will be: (1) a search for free- free emission from ionized shells surrounding these clusters 
at mm wavelengths, permitting a clear discrimination between thermal and non-thermal 
emission; (2) a search for near-IR emission with colors signifying the stellar content of the 
SSC. Such continuum emission would not be espected from the remnant of an isolated SN 
since any associated cluster would produce a SN wind driven supershell rather than an 
isolated remnant; (3) near IR measurement of the ratios of [Fell] to recombination line 
emission; and (4) measurement of shell expansion velocities. These tests are feasible now 
with deep searches and radio observations at high angular resolution made over a long time 
baseline. 

The WDB picture needs to be harmonized with the larger picture of M82. For example, 
why aren't young SSCs with masses outside the range indicated in Figures 2 and 3 also 
associated with non-thermal sources? In the WDB model, should we see an association 
between optically visible SSCs and CNRs? If the CNRs are not SNRs then where are the 
latter objects? We consider these issues in turn. 

Naturally, SSCs with masses lower than the observed range would not have winds pow- 
erful enough to yield sufficient relativistic particle energy. Clusters with higher masses could 
emit synchrotron radiation, but such emission could be suppressed by inverse Compton cool- 
ing if the ambient gas density is high enough to confine the bubble for a sufficiently long 
period, as shown in section 5.1. 

Regarding the potential association between CNRs and v isible SSCs, such identification 



is unlikely, since these clusters are found to have ages > 10 6 yr (jMelo et al.ll2005l ). and are thus 
too old to exhibit nonthermal emitting shells. The very young clusters associated with this 
CNR phenomenon would very likely be obscured by the surrounding dust in the molecular 
clouds where they are formed. In addition, the current lack of astrometric precision for 
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optically visible clusters precludes a serious search for identification, especially for the faint 
members with lower masses (near 1O 4 M ). 

Concerning the SNRs in M82, SNe may originate predominantly in SSCs if such clusters 
are the most common sites of massive star formation. Individual SNRs may then be more rare 
than previously assumed, since the aggregate effect of SNe ejecta in clusters is to produce SN 
winds even more powerful than the earlier stellar winds. The SN winds would re-energize the 
pressure confined shells produced by the earlier winds to produce the supershells evident in 
M 82 by free-fre e and HI absorption of the non-thermal background and by CO emi ssion from 
molecular gas Jwills et al.lll997l : IWeiB et al.l ll999l : IPedlar. Muxlow fc Willsl l2003h . Though 
such shells may also produce non-thermal emission, their size (tens to hundreds of pc) would 
preclude easy identification because of the general confusion by the complex non-thermal 
background. It is undoubtedly such shells, particularly from massive clusters, which break 
out of the disk and contribute to energizing the nuclear outflows in M 82 and other starburst 
galaxies. Note that if many of the CNRs are not SNRs, then the SNe rate based on CNRs 
alone would be overestimated. However, this deficit may be more than compensated by the 
complementary rates inferred from the interpretation that supershells are the product of 
multiple SNe in SSCs. 

The compact non-thermal sources in M82 and other galaxies may thus be the unique 
signatures of the winds of low mass very young clusters embedded in the ISM, except for 
some individual cases where their time variability and/or measured expansion rate identify 
them uniquely as SNRs. 
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A. Appendix: [Fell] Emission from Cluster Wind Driven Bubbles 



There are two potential sources of [Fell] emission from a WDB. The first is collisionally 
excited gas in the radiative shock formed as the interstellar medium is swept up into a 
shell by the expanding bub ble. The second is the part of this shell photoionzed by the 
cluster. iMouri et al.l (120001 ) has modeled [Fell] 1.257/im emission from both shocked and 
photoionized gas, and we base our estimate on this work. The results are applicable as 
well to emission at 1.644 /im, relevant to the discussion in §8.1, since both transitions have a 
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common upper state and have very similar transition probabilites (e.g. iNussbaumer fc Storey 
1988h . 



(a) Emission from shocked gas in the shell 

The [Fell] luminosity L s from the shocked gas may be parameterized as a fraction 
?7 S = L s /Efc of the flow rate of kinetic energy of ambient gas into the shock, which is 
al so the rate o f energ y radiated by the shell. The parameter 7] s is gi ven in Figure 11(b) 
of iMouri et all (bood ). The quantity E fc is given by E fc = (27/77) L w jMac LowfeMcCrav 



19881 ). Combining this with equation (3b) gives 



L s = 2.0 x 10 37 r7 s M 4 ergs~ 1 
(b) Emission from photoionized gas in the shell 



(Al) 



The [Fell] luminosity L p is estimated here for an ionization bounded shell, and may be 
written as L p = (L p /Pa/3)(Pa/3/Ly)Ly, where Pa/3 and Ly are respectively the Paschen f3 
and Lyman continuum luminosities. The ratio r] p = (L p /Pa/5) is dependent on the assumed 
metal abundances, the shell gas density, the spectru m of the sour c e of i onization, and the 
ionization parameter. It is given by Figure 8(a) of IMouri et al.l (120001 ) for a limited set 
of co nditions. The ra tio (Pa/3/Ly) may be computed from standard case B recombination 
(e.g. lOsterbrockl Il989l ) using an electron temperature of 10 4 K. Combining this result with 



equation (3e) then gives 



L p = 3.3 x 10 37 r7 p M 4 ergs~ 1 



(A2) 
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